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EFFECTS OF EDDY VISCOSITY ON TIME CORRELATIONS 
IN LARGE EDDY SIMULATION 


GUOWEI HE* , R. RUBINSTEIN*, AND LIAN-PING WANG* 

Abstract. Subgrid-scale (SGS) models for large eddy simulation (LES) have generally been evaluated by 
their ability to predict single-time statistics of turbulent flows such as kinetic energy and Reynolds stresses. 
Recent applications of large eddy simulation to the evaluation of sound sources in turbulent flows, a problem 
in which time correlations determine the frequency distribution of acoustic radiation, suggest that subgrid 
models should also be evaluated by their ability to predict time correlations in turbulent flows. 

This paper compares the two-point, two-time Eulerian velocity correlation evaluated from direct numerical 
simulation (DNS) with that evaluated from LES, using a spectral eddy viscosity, for isotropic homogeneous 
turbulence. It is found that the LES fields are too coherent, in the sense that their time correlations decay 
more slowly than the corresponding time correlations in the DNS fields. This observation is confirmed by 
theoretical estimates of time correlations using the Taylor expansion technique. The reason for the slower 
decay is that the eddy viscosity does not include the random backscatter, which decorrelates fluid motion at 
large scales. An effective eddy viscosity associated with time correlations is formulated, to which the eddy 
viscosity associated with energy transfer is a leading order approximation. 

Key words, eddy viscosity, time correlation, large eddy simulation 

Subject classification. Fluid Mechanics 

1. Introduction. In large eddy simulation of turbulent flows, large-scale velocity fields are computed 
directly from the filtered Navier-Stokes (N-S) equation, while the effects of small-scale velocity fields on 
large-scale velocity fields are modeled using a SGS model. The SGS models are constructed to represent 
energy drain from large scales to small scales, and local inverse energy transfer, i.e., energy backscatter (from 
small scales to large scales). It is desirable that the SGS models can produce a LES field which correctly 
predicts the large-scale statistics of the N-S field. A direct evaluation [1, 2, 3] of the SGS models is to compare 
statistics of the LES fields with those of the filtered N-S fields. There has been extensive work comparing 
single-time statistics of the LES and DNS fields. For recent reviews, see [4, 5]. However, there is little work 
concerning two-time statistics. This paper investigates the effects of the eddy viscosity SGS models on the 
time correlations of the LES fields. The research is motivated by use of LES in the aeroacoustics [4, 6], where 
sound radiation is dependent on Eulerian time correlations. 

Time correlations are among the simplest statistical properties of turbulent flow. We will consider 
Eulerian time correlations in isotropic homogeneous turbulence, whose Fourier transformation is expressed 
as 


(i.i) 


C(k , r) = (u(k, t + t) ■ u(k, t)). 
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where u(k, t) is a Fourier mode of the velocity field at the wavenumber vector k with k = |k|. The bracket 
“< )” means ensemble averaging. The Eulerian time correlation measures the temporal changes seen by 
an observer in a fixed reference frame with zero mean velocity. For a possible equilibrium range at high 
frequencies, such an observer sees small scale eddies being swept past him by larger scale eddies. The 
correlation coefficient [8, 9] between the LES and DNS fields has been used to evaluate the ability of LES 
to predict flow evolution. It was demonstrated that the correlation coefficient will decay to zero after a few 
eddy turnover times. This implies the LES and DNS fields will be completely decorrelated. However, this 
result is not discouraging. In fact, two DNS fields which are initially identical at large scales but slight 
different at small scales will also become completely decorrelated, although their uncoupled statistics still 
remain equal. Even independent random fields can be statistically identical. Therefore, the relevant question 
to ask is if the LES fields can reproduce the statistics of the DNS fields, whether the LES and DNS fields are 
correlated or not. In this paper, the validity of the LES is evaluated by its ability to reproduce the Eulerian 
time correlations of the DNS field. 

Current SGS models can be categorized into several groups: eddy viscosity, stochastic, similarity, as- 
sumed SGS velocity models [5] and multiscale model [10]. The most extensively used SGS model is based 
on the eddy-viscosity assumption: the effects of subgrid scale motion is supposed to be similar to that of 
molecular dissipation, but with an eddy viscosity which depends on both wavenumber and the cutoff. Many 
eddy viscosity models have been developed including the Smagorinsky model, the kinetic-energy model and 
the dynamical model. In isotropic homogeneous turbulence, the eddy viscosity can be formulated in terms 
of the energy balance equation: 

(i.2) f^ + 2i/* 2 Jjs(fc,o = r(M), 

where E{k,t) is the kinetic energy spectrum, v the kinematic viscosity and T(k,t) the energy transfer. One 
introduces a cutoff wavenumber k c such that only scales larger than 1 jk c will be resolved explicitly. The 
energy transfer can be separated into two parts 


(1.3) T(M) = T<(k\k e ,t) + T>(k\k c ,t), 

where T < (k\k c , t) is the part of the total energy transfer induced by triad interactions among waves p, q and 
k satisfying p-f q = k such that p and q are less than k c , and T > (/r|fc c , t) the part of the total energy transfer 
induced by triad interaction among the wavenumber p, q and k satisfying p-fq — k such that p or q is larger 
than k c . Hereinafter, the superscript < implies the part of any quantity involving only wavenumbers less 
than the cutoff k C7 while the superscript > denotes the remainder, in which at least one of the wavenumbers 
is larger than the cutoff k c . The eddy viscosity is defined by [11]: 


(1.4) 

so that finally 


v t {h\k c ,t) - - 


T>(k\k c ,t) 
2 k 2 E{k,t) 


(1.5) 


+ 2{v + v t )k 2 


E(k,t) = T < (k\k c ,t). 
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In the formulation of the papers [12, 13] 


(1.6) 


Vt{k\k c ,t) = vf(k\k c ) 


I E(k c ,t) 


where E(k c ,t) is the energy spectrum at the cutoff wavenumber k c and vf (fc|/c c ) is 


(1.7) Vt{k\kc) = 0.267 + 9.21 exp(-3.03Jfc c /fc). 

This eddy viscosity model has been extensively studied by DNS and experiments. These studies show 
that it does model the energy dissipation at small scales but misses the energy backscatter from small scales 
to large scales. We will show that the lack of the random backscatter causes more coherence, in the sense 
that the time correlations of the LES fields decay more slowly than those of the DNS fields. We will focus 
on the eddy viscosity model (1.4). The results may be extended to other eddy- viscosity SGS models. In 
Section 2, we calculate the Eulerian time correlations for the DNS and LES fields, respectively. We will show 
that the LES fields predict larger time correlations. In Section 3, the reasons for the larger predictions are 
explained. We use the Taylor expansion technique to estimate the time correlations respectively for the N-S 
and LES fields. The effective eddy viscosity for time correlation is derived. Discussion and conclusions are 
given in Section 4. 

2. Numerical calculation of time correlations. We carry out DNS and LES of isotropic homoge- 
neous turbulence. The three-dimensional N-S equations are solved numerically in a periodic box of side 27r 
by the standard pseudospectral algorithm. In DNS, the flow cube is discretized uniformly into A r = 128 grid 
points, which defines the wavenumber components in Fourier space as kj == ±0, 1, ■ • ■ , N/2 — 1 for j = 1,2, 3. 
The aliasing errors are removed by the two- thirds truncation method. Consequently, the effective wavenum- 
ber space is by \k\ < N/ 3. The initial condition is set to be isotropic and homogeneous with the energy 
spectrum E(fc, 0) oc k 4 /&q exp(— 2(fc/fco) 2 ), where ko is the wavenumber at which the maximum of the energy 
spectrum occurs. By steady random forcing on the first two wavenumber shell k < 2, stationary conditions 
are achieved after some passage of time. In this study, the relevant parameters are: viscosity v = 0.002, time 
step size dt — 0.005, Taylor microscale wave number R\ = 40, CFL number 0.46 and velocity derivative 
skewness —0.34. In LES, N = 64. We use the eddy viscosity model (1.6) and the sharp cutoff filter. All 
other parameters in the LES are kept to be the same as the DNS. The spectral codes for the DNS and LES 
are developed and implemented on the ICASE Beowulf 96 CPU cluster computer, Coral. 

In Fig. 2.1, we plot the time correlations of the LES fields against the time delay. Fig. 2.1 clearly 
demonstrates the spatial scale dependence of the time correlations: the larger wavenumbers remain coherent 
for a longer time while the smaller wavenumbers are rapidly decorrelated. Fig. 2.2 and Fig. 2.3 show the 
same data against time delay normalized by eddy turnover time r e = [/?a/A:E(^)] _ 1 and sweeping time 
r s = ( uk)~ l respectively, where u is the rms velocity. The sweeping time is more suitable for collapsing the 
time correlation. 

Fig. 2.4 shows the time correlations of the DNS and LES fields for different wavenumbers, k = 5, 10, 15, 20, 
spanning a range from the scale on the order of the integral scale to the scale on the upper end of resolved 
scale range. The times are normalized by the sweeping time at the largest wavenumber k — 1. It is evident 
that the time correlations of the LES fields decay more slowly than those of the DNS fields. 
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Fig. 2.4. Time correlation C{k,r) versus the time lag t/t s . The solid line for the LES field and the dash line for the 
DNS field: k = 5,10,15,20. 


3. The effects of eddy viscosity on time correlation. The numerical observations in the last 
section can be supported by theoretical estimates of time correlations. We will use the Taylor expansion 
technique to estimate time correlations for both the N-S and LES velocity fields. The Taylor expansions have 
been developed to construct approximations to statistics in turbulence. Kaneda’s group [14, 15] has done 
extensive work on the Taylor expansion approach. Its central idea is to expand the time correlation into a 
power series of time lag r. The coefficients in the Taylor series are determined by the governing equations. 
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The governing equations for the N-S fields it Q , a = 1, 2,3, in spectrum space are 


(3.1) ^ + vk-j « a (k, t) = M a fo (p + q = k)u J u(p)« 7 (q), 
where 

M tt/ j 7 (p + q = k) = Af a/ , 7 (k) ^ > 

P f q=k 

M«/J 7 (k) = Yi {kj)P a -y{ k) + k 1 P n0 (k)) , 

(3.2) F> a ^(k) = <5 a> J - kakp/k 2 , 

and Sp+ q =k denotes a summation for all p and q such that p + q = k. 

The governing equations for the LES fields u a , a = 1,2,3, in spectrum space are 

(3.3) {§i + [ l ' + l ' t (*lk)l* 2 ) «« (k, t) = (p + q | = k)w.j(p)w 7 (q) . 

To facilitate comparison, we simply assume that the N-S and LES velocity fields u a (k, t) are identical at the 
time t. The velocity fields u Q (k, t + r) at the later time t + r are obtained from the equation. 

The Taylor expansion of the time correlation (1.1) has the following form 


C(k, r) = C 0 (k) + Ci ( k)r + C 2 (k)r 2 / 2 + ■ ■ ■ , 

■ (3.4) C n (k) = (^^-n(-k,t}). 

The time derivatives of the velocity fields can be formally calculated from the N-S equation. For a negligible 
viscosity v } we obtain 


C 0 (k) = <n Q (k,£)ua(-M)), 

c.(*) = |(r< + r>), 

(3.5) C 2 (k) = 2M at jy(p + q = k)A/^ cr (m + n = p)(u M (m,t)u CT (n,«)u 7 (q,t)u a (-M)). 

Similarly, the Taylor coefficients of time correlations in the LES velocity fields can be obtained 


C 0 (k) - (u Q (k, t)u Q (— k, t)), 

Ci{k) = 1 (Tl r Jl/|(*|*e)k a <«a(k > 0«a(-k,t))) , 

C 2 {k) = 2M < - Mai}-, (p + q = k)/y,(fc!C)(p 2 + q 2 = k 2 ){u 0 (pj)u-,(q,t)u a (-k, <)> 

(3.6) +[is t (k\k c )k 2 } 2 (u a (k,t)u a (-k,t)}. 

Comparing the equations (3.5) with the equations (3.6), we find: (1) Co(k) = Cq (At); and (2) the eddy 
viscosity (1.4) implies Ci(k) = C\ (k). In other words, the eddy viscosity (1.4) can be obtained solving 
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the equation C\ (k) = C\(k). Recall that Ci(k) = 0; (3) Cz(k) 7^ (?2(fc) for the eddy viscosity (1.4). This 
implies that the time microscale of the LES field is not equal to that in the DNS field. Therefore, the time 
correlations of the LES fields are different from those of the N-S fields. The numerical calculation indicates 
that the former is larger than the latter. 

The results from numerical calculation and theoretical estimations can be understood by the following 
physical arguments: The contributions of small scales to large scales can be described as energy dissipation 
and random backscatter. The eddy viscosity correctly models the drain of energy from large scales to small 
scales but misses the random backscatter from small scales to large scales. This leads to a more coherent LES 
field. Therefore, the LES field evolves in a more correlated fashion, in the sense that their time correlations 
decay more slowly. 

The eddy viscosity is introduced as a surrogate for energy drain from resolved scales to the subgrid scales. 
Its expression (1.4) is then found from the energy balance equation. Following this idea, the effective eddy 
viscosity for time correlation can be found from its governing equation as a surrogate of the contribution 
from subgrid scales to resolved scales. 

The governing equation for time correlation may be written as 

(Jr + "* 2 ) + r ) u a(“M)> 

(3.7) = AWp + q = k) ((u d (p, t + r)tt-y(q, t + r),u Q (-k,£ )) - c.c) 

where “c.c” represents the complex conjugate of the triple moment. We separate the nonlinear interaction 
on the right hand of the equation (3.7) into two parts: for p and q < k and for p or q < k such 

that — M<3 -f . Then, the effective eddy viscosity which describes the effects of the subgrid 

scales on the time correlations of the LES fields may be defined by 


(3.8) 


v T {k\k c ) 




2k 2 (u a (k, t + r)u Q (-k,£))* 


We can also derive a Taylor expansion of the effective eddy viscosity v T {k\k c ) in power of time lag r, leading 
to 


(3.9) v T {k\k c ) = v t (k\k c ) + O(r). 

It indicates that the normal eddy viscosity (1.4) is the leading order approximation to the effective 
eddy viscosity (3.8). For small time lag, the normal eddy viscosity produces a satisfactory approximation. 
However, for large time lag, the corrections from higher order terms must be taken into account. 

4. Discussion and conclusion. The normal eddy viscosity models the energy dissipation but ignores 
the random backscatter. It produces a more coherent LES field. It is shown by numerical calculation and 
theoretical estimation that the time correlations for the LES field decay more slower than those for the 
N-S field. The differences in time correlations between the LES and DNS fields decrease with increasing 
wavenumber. The effects of the eddy viscosity on time correlation have to be considered while applying LES 
to sound radiation. 

We have shown that the normal eddy viscosity for the energy balance equation is a leading order 
approximation to the effective eddy viscosity for the time correlation equation. Therefore, it is difficult to 
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construct a new surrogate of eddy viscosity, which could exactly satisfy both energy balance equation and time 
correlation equation. An effective eddy viscosity that could simultaneously produce a good approximation to 
the terms involving small scales in the energy balance and time correlation equations will be helpful. Noting 
that the information given by the eddy viscosity are only statistically averaged one, we need to account 
for the stochastical effects on the particular realization of the LES fields. The stochastical effects could be 
introduced by random forcing. Therefore, the decorrelation of the LES field may be resorted by including 
random force. 
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